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Abstract 

We discuss the diffractive deep inelastic scattering at very small x and derive the properties 
of the diffractive dissociation of virtual photons in the triple-pomeron regime. We demonstrate 
that the photon-pomeron interactions can be described by the partonic structure function, which 
satisfies the QCD evolution equations, and identify the valence and sea (anti) quark and the va- 
lence gluon structure functions of the pomeron. The gluon structure function of the pomeron 
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can be described by the constituent gluon wave function of the pomeron. We derive the leading 
unitarization correction to the rising structure functions at small x and conclude that the unita- 
rized structure function satisfies the linear evolution equations. This resul holds even when the 
multipomeron exchanges are included. 
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1 Introduction. 



The pomeron (IP) remains one of the most misterious objects in the high energy physics. Apart 
from the elastic scattering, the exchange by pomerons describes (Fig. 1) the diffraction dissoci- 
ation of the projectile, which can be viewed as the projectile-pomeron interaction (Fig. Ic) [1]. 
In the diffractive leptoproduction at x = Q"^ /{Q"^ + W^) <ti 1 one can think of the deep inelastic 
scattering (DIS) on the pomeron emitted by the target nucleon [2-8]. (Here W is the total 
energy in the photon-proton center of mass, = 2pq — Q^, where p and q are the 4-momenta 
of the proton and photon, and = —q^ is the virtuality of the photon.) If the diffraction 
dissociation is dominated by the single pomeron exchange, which is a very strong assumption, 
and if the pomeron exchange can be treated as the factorizing particle exchnage, which also 
is a very strong assumption, then one can introduce the operational definition of the (virtual) 
photon-pomeron cross section (Ttot(7*IP, Q^, M^) and of the structure function of the pomeron 
F2^\x,Q'^) in terms of the differential cross section da d / dtdM'^ of the forward diffraction dis- 
sociation of virtual photons 7* + p — X + p ( we follow the Regge-theory convention [1] with 
the substitution + Q"^ which is natural for the DIS): 



a,„,(7*IP, M^) = ^^(M^ + Q^^ +P^X+p) 



atotipp) dtdM^ 
and 



t=0 



FPix, Q') = -^atoth*W, M') (2) 
with the corresponding Bjorken variable 

There has already been much work on the parton model phenomenology of the pomeron 
[2-10], but the definitive proof that the so-defined structure function of the pomeron satisfies the 
conventional Gribov-Lipatov-Dokshitzer-Altarelli-Parisi (GLDAP) QCD evolution equations [11- 
13] is as yet lacking. The definition (P for atotil*'^) does implicitly assume that the pomeron 
has the intercept q;ip(0) = 1, i.e., the high energy cross sections are constant and the mass 
spectrum of excitation of large masses M has the 1/M^ behavior. 



If the factorization relations are valid, then A^jp is expected to be a universal dimensional con- 
stant independent of the projectile a. However, in QCD there are no a priori reasons for the 
factorization relations to hold, and there are indications to the contrary [7,8,14,15]. Further- 
more, the factorization could strongly be violated by the absorption (unitarity) effects from 
the multiple pomeron exchanges in Figs, lb, Id [6]. These multiple pomeron exchanges cast 
shadow on the reinterpretation of the diffraction dissociation in terms of the photon-pomeron 
interaction. Experimentally q;ip(0) > 1 [16,17], the total cross sections are rising and the mass 
spectrum of the diffraction dissociation of protons exhibits slight deviations from the 1/M^ law 
[18]. Furthermore, the cross section of the diffraction dissociation of virtual photons was shown 
to be infrared sensitive [7-9,19]. The quantity related to the diffraction dissociation cross section 
- the unitarization (shadowing, absorption) correction to the structure functions at small x - is 
also infrared sensitive [19-22] . Therefore, the possibility of introduction of the well-defined, and 
well-behaved in the sense of the QCD evolution, structure function of the pomeron and the issue 
of the infrared sensitivity of the diffraction dissociation cross section deserve further study. 

In this paper we address the above problems in the framework the light-cone s-channel 
approach to the diffractive deep inelastic scattering at small x = Q'^/(W^ + Q"^) initiated in 
our previous publications [7,23,24] (for the related early work on the s-channel approach to the 
light-cone QED see Bjorken, Kogut and Soper [25]). Our strategy is to compute the high energy 
behavior of the total (virtual) photoabsorption cross section and of the diffraction dissociation 
cross section. We treat the diffractive 7*p scattering in terms of absorption of the light-cone 
partonic Fock components of the virtual photon on the target proton or nucleus. One can do so 
since at a; <^ 1 the photon transforms into these partonic Fock components at large distances 

~ ^— > (5) 

upstream the target nucleon (nucleus). As an illustration of the principal points of the light- 
cone formalism [24,7] consider interactions of the qq Fock state of the photon. Because of the 
condition (|^) the transverse size p of the qq pair and the partition z and {1 — z) of the (light- 
cone) momentum of the photon between the quark and antiquark can be considered frozen in the 
scattering process. Therefore, one can introduce the spatial wave function of the light-cone qq 
Fock states ^^*{f>, z) and the dipole cross section cr(p) such that the total photoabsorption cross 
section aT,L{l*N,x,Q^) for the (T) transverse and (L) longitudinal photons and the forward 
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diffraction dissociation cross section can be calculated as the conventional quantum-mechanical 
expectation value [24,7] 

aT,L{l*N,x,Q^)= £dz J dVl^T,L(^,p)|V(p) , (6) 

= Tdz [d'p\<ilTA^,p)\Mpr ■ (7) 

iOTT JO J 

We emphasize that the factorization of the integrands in Eqs. is exact, and corresponds to 
the exact diagonalization of the scattering matrix in the (p, z)-representation. 

Interactions of the qq Fock state of the photon give the driving terms of the structure function 
at small x and of the diffraction dissociation cross section. Specifically, Eq. @ yields the 
photoabsorption cross section and the proton structure function F2^\x, Q"^) which are constant 
vs. X. Eq. yields the mass spectrum of the diffractively produced states daoi^* ~^ X)/dM^ oc 
+ M"^)"^ and can be associated with the 'valence' qq component of the pomeron [7,23,24]. 
Both the rise of F2^\x,Q'^) towards small x and the triple-pomeron component of the mass 
spectrum in the diffraction dissociation of photons, which can be associated with the 'sea' qq 
pairs and gluons in the pomeron, are generated by interactions of the higher, qqgi....gn Fock 
states of the photon. 

The subject of this paper is a generalization of the light-cone s-channel approach [24,7] to 
interactions of the higher Fock states of the photon. Our major findings can be summarized as 
follows: The diffraction dissociation cross section can indeed be factorized (Eq. (PP|)) into the 
flux of pomerons in the proton f-sp{y)/y, where y = (M^ + Q^)/ {W'^ + Q'^) is a fraction of proton's 
(light-cone) momentum carried by the emitted pomeron, and the well-defined structure function 
of the pomeron F2^\x/y,Q'^) , which satisfies the conventional QCD evolution with (the 
definition ([1|) corresponds to the convention /ip(y) = 1). This structure function and the flux of 
pomerons describe how the naive oc mass spectrum (^ is modified by the rising hadronic 

cross sections and by the QCD evolution effects. The factorization ( P^ ) bears certain semblance 
to the usual Regge-theory factorization despite the fact that the Regge-theory factorization 
relations do not hold in deep inelastic scattering at small x and in our analysis we never assume 
nor use the Regge-theory factorization. We speak of the triple-pomeron regime just to pay a 
tribute to the fact that we mostly consider <^ <^ W"^, which in the Regge-theory would 
have been the triple-pomeron domain. The infrared sensitivity of the diffraction dissociation 



dt 



t=o 



dM' 



2 daoil* 



X) 



dtdM^ 
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cross section can be reabsorbed into the initial momentum distribution of the 'constituent' 
quarks, antiquarks and gluons in the pomeron, in close similarity to the conventional QCD 
evolution analysis of the proton structure function. We demonstrate that besides the 'valence' 
and 'sea' quark-antiquark components derived in [7] (for discussion of the valence qq component 
of the pomeron see also [5]), the pomeron has the 'valence' gluon component, and present 
the explicit derivation of the constituent gluon wave function of the pomeron. The absolute 
normalization of the pomeron structure function can be related to the triple-pomeron coupling 
^3ip(Q^) Eq. ( |62D , which gives the driving term of the diffraction dissociation mass spectrum (^. 
We shall confirm the earlier suggestion [26] that despite being dimensionfull, Asip{Q'^) oc GeV~^, 
this triple-pomeron coupling only has weak dependence on Q^, although if were the only 

relevant scale in the DIS, then naively one could have expected ^3ip(<5^) oc l/Q^. The fact, 
that one can (approximately) relate the properties of the diffraction dissociation in the DIS and 
in the real photoabsorption, which here is proven rather than assumed, is interesting by itself. 
In [7] we gave a phenomenological estimate of the sea structure function of the pomeron and of 
the diffraction dissociation rate in terms of the triple-pomeron coupling A3]p(0) borrowed from 
the real photoproduction data [27]. The resulting predictions are in a good agreement with the 
first data on the diffraction dissociation of photons in DIS at HERA obtained recently by the 
ZEUS collaboration [28]. 

The diffraction dissociation of virtual photons and the unitarity (shadowing, absorption) 
corrections to the rising structure functions at small x are two closely related phenomena. The 
GLDAP evolution equations predict [12] very rapid rise of the structure functions towards large 
which conflicts the s-channel unitarity [19]. The s-channel unitarization of the virtual 
photoabsorption cross section introduces the multiple-pomeron exchanges, which could lead to 
the departure of the x- and Q^-dependence of structure function from predictions of the GLDAP 
evolution. We find that the unitarity corrections are large and persist at all Q^. The unitarity 
correction is a nonlinear functional of the parton density and violates the conventional linear 
relationship between the photoabsorption cross section and the parton density. Our principle 
conclusion is that, nonetheless, the modified evolution equations retain their linear GLDAP 
form as distinct from the nonlinear equation suggested in [19] and discussed extensively in 
the literature during the past decade (for a recent review with many references see [22]). The 
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diagonalization of the scattering matrix in the (p, ^)-representation greatly simpHfies a discussion 
of the unitarization correction, as it allows to unambiguously identify the s-channel partial waves 
which must satisfy the unitarity bound. 

This paper is organized as follows. In section 2 we review the light-cone s-channel approach to 
deep inelastic scattering starting with the diffractive interactions of the two-body qq Fock state 
of the photon [24,7]. In section 3 we study interactions of the 3-body qqg Fock state of the photon 
and derive the driving term of the triple-pomeron mass spectrum in the diffraction dissociation of 
virtual photons and the driving term of the sea structure function of the pomeron. In section 4 we 
apply our formalism to derivation of the rising structure function in the double- logarithmic limit. 
The subject of section 5 is the triple-pomeron regime of the diffraction dissociation of photons to 
higher orders in the perturbative QCD. Here we derive the structure function and the constituent 
gluon wave function of the pomeron. The latter absorbs the infared-regularization dependence 
of the diffraction dissociation cross section. This completes a derivation of the valence and sea 
(anti) quark and the gluon structure functions of the pomeron, which are to be used as an input 
of the GLDAP evolution of the pomeron structure function. The factorization properties of the 
QCD pomeron and the flux of pomerons in the proton are discussed in section 6. In section 
7 we discuss the unitarization of the rising structure functions of the proton at small x. We 
demonstrate that the unitarized structure functions of the proton still satisfy the hnear GLDAP 
evolution equations, as distinct from the nonlinear Gribov- Levin- Ryskin (GLR) equations [19]. 
Remarkably, this conclusion holds even when the multiple pomeron exchanges are included. In 
this section we also present a brief phenomenology of the shadowing correction to the proton 
structure function and comment on the treacherous path to the interpretation of the shadowing 
in terms of the fusion of partons. In section 8 we comment on the relation between our scenario 
for deep inelastic scattering at small x and the Lipatov's work on the pomeron in QCD [29]. 
Our main results are summarized in section 9. 

This paper is mostly devoted to the derivation of the formalism, the numerical results for 
the unitarization of the proton structure function at small x [30], for the (anti)quark and gluon 
distributions in the pomeron and the structure function of the pomeron [31], for the fusion of 
the nuclear partons and the nuclear shadowing [32] will be presented elsewhere. 
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2 Deep inelastic scattering in terms of the Fock states of the photon 
and the dipole cross section 



2.1 The qq Fock states of the photon and sea quarks in the proton 

We are interested in DIS at x ^ 1 where the structure functions are dominated by the scattering 
of photons on the sea quarks. The driving term of the sea is given by the perturbative QCD 
diagrams shown in Fig. 2. The same diagrams can be viewed as the scattering on the target proton 
of the qq Fock states of the photon. The principle finding of Ref. [24] is that the corresponding 
contribution to the photoabsorption cross section can be cast into the quantum-mechanical form 
d^). The wave functions of the qq fiuctuations of the photon were derived in [24] and read (for 
the related discussion in the framework of QED see also Bjorken et al. [25]) 

I^t(.,p)P = ^,EZj{[z' + (1 - zryK^iepr + rajMepf} , (8) 

\^Liz,p)\' = ^E4^? Q'^'i^ - ^fM^Pr , (9) 

where Ki,(x) are the modified Bessel functions and 



e'^ = z{l - z)Q'^ + mj . (10) 

In Eqs. (pD-(p!OD is the quark mass and z is the Sudakov variable, i.e. the fraction of photon's 
light-cone momentum g_ carried by one of the quarks of the pair (0 < 2; < 1). In the diagrams of 
Fig. 2 the color-singlet qq interacts with the target nucleon via the Low-Nussinov [33] two-gluon 
exchange, which is the driving term of the QCD pomeron. The interaction cross section for the 
color dipole of size p is given by [24] 

16 , , f d'^kV{k)[l-exp{ikp}] --^ 
'Tip) = -jasip) J la.(P) , (11) 

In ( pT| ) k is the transverse momentum of the exchanged gluons, the longitudinal momentum of 
the exchanged gluons is ~ m^x and is negligible at small x, pc ~ V-^c is some kind of effective 
mass of gluons introduced so that color forces do not propagate beyond the confinement radius 
Rc ~ Rn- The gluon-gluon-nucleon vertex function V{k) is related to the two-quark form factor 
of the nucleon (j2(^i, ^2) = {N\ exp i{ki ■ ri + ^2 ■ ^2) |^) by 

vik) = 1 - G2{k, -k)^i- :f,^{?,P) , (12) 



where Tchi^f) is the charge form factor of the proton. ^^(A;^) is the running QCD couphng 
which we shall use both in the momentum, and the coordinate representations: 

An An 

^s{k^) = — n2/\2 V' ^s{p) = — 777T2 ■ (13) 

/3o log(fc7A^c^) polog{l/ A'qcdP^) 

2.2 Universality of the dipole cross section and infrared regularization 

The salient feature of the dipole cross section (|ll|) is its universality property [7,8]: cr(p) de- 
pends only on the size p of the qq color dipole. The dependence on and the quark flavor is 
concentrated in the wave functions @ and (j^). The fundamental role of the color gauge invari- 
ance must be emphasized. By virtue of the color gauge invariance gluons with the wavelength 
\ = 1/k > Rj^ decouple from the color-singlet nucleon. This decoupling is taken care of by the 
vertex function V{k), which vanishes as — > 0, and the size of the nucleon emerges as a natural 
infrared regiularization: the dipole cross section (0) is infrared-finite even if fie = 0. Similarly, 
the factor [l — exp{ikp)] takes care of the decoupling of gluons with X > p from the colour-singlet 
qq Fock state. As a result, at small p the cross section a{p) is perturbatively calculable with its 
absolute normalization. For the nucleon target 



, , 47r 2 , , f^p' dkW{k) .j^^. 16n' , , , , 
^(P) ~ -irP as{p) / -jr- ^as{k') = ——p' as{p) log 



1 



where 



Lip) = log 



CNp'as{p)L{p),{lA) 



(15) 



as{p)_ 

is the large parameter of the so-called Leading-Log Approximation (LLA) [11]. 

Another universal feature of a{p) is its saturation at p > Rn, Rc because of the confinement 
[24]. This is a strong-coupling regime and in the regime of saturation cr(p) depends on the infrared 
regularizations. (Following [34], we assume the freezing strong coupling as{p) = as{Rf) ~ 1 at 
p > Rf.) One natural infared regularization - the size of the target proton - enters via the vertex 
function V{k). The other two infrared regularizations are the effective confinement radius Rc 
which enters via the effective gluon mass pc and the freezing point Rf of the strong coupling 
as{Rf) ~ 1- Evidently, by virtue of Eq. (^ the dependence on these infrared regularizations 
propagates into the proton structure function. Here we would like to emphasize that such an 
infrared sensitivity of F2^\x, Q"^) is old news: in the conventional QCD-improved parton model 
this dependence is hidden in the parameterization of the input parton distributions at the low 
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factorization scale Q"^. In our light-cone s-channel approach we rather calculate the structure 
function in terms of the dipole cross section, reducing drastically the number of the infrared 
parameters. Furthermore, the crude test of the large-p behavior of the dipole cross section (j{p) 
is provided by the hadronic cross sections. For instance, the pion-nucleon total cross section can 
be evaluated as 

Ototi-r^N) = I'dzJ d'p |^,(z,p)|V(p) , (16) 

which roughly reproduces the observed value of atot{.T^N) at moderate energies [7,8,14,15,30]. 
Once the constraint ( |16|) has been enforced, the predictions for DIS become to a large extent 
parameter-free. Such a minimal-regularization approach leads to a viable description of the 
absolute value of F2^\x, Q^) [23,24,30,35], of the longitudinal structure function Fl{x, Q"^) [36] 
and of the nuclear shadowing in DIS [23,24,35], of the gluon distribution in the proton [37], of 
the excitation of charm in the muon and neutrino scattering [38] and roughly reproduces the 
total cross section of the real photoabsorption [7,8]. This minimal-regularization approach is 
not imperative, though, and all the results to be presented below could easily be reformulated 
in terms of more familiar parameterizations of the input parton distributions in the proton and 
pomeron at the expense of losing the predictive power. 

2.3 Connection with the QCD evolution equations 

At small the transverse size of the qq Fock states of the photon shall be given by the Compton 
wavelength of the quark l/rrif. For the heavy flavors (charm,...) this is the perturbatively small 
size. For the light flavors it is natural to ask that also the quarks do not propagate beyond the 
confinement radius. With the natural choice rriu^d ^ fJ'n ^ ^/Rc Eq- (§|) with the wave function 
(H) reproduces the real photoabsorption cross section [7,8]. At larger ^ mj, Pq, R]^^ the 
conventional QCD-improved parton model description is recovered. Indeed, let us calculate the 
dependence of the cross section (|]). At large the leading contribution to aT{'y*N) comes 
from the Ki{epY term in Eq. (^. Making use of the properties of the modified Bessel functions, 
after the 2;-integration one can write 

aT{i*N.x,Q') = fdzf dVl*T(^,p)|V(p) oc f'^^ ^aip) 

Jo J J1/Q2 p^ 

-^]^^2^L{Qr. (17) 



^7yl , o —^s{p) log 
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We find the scaling cross section oc l/Q^ times the LLA scahng violation factor, with one power 
of L{Q'^) = logll/asiQ"^)] per QCD loop, which is a starting point of the derivation [11-13] 
of the QCD evolution equations. Notice, that the factor l/Q"^ in Eq. (|l^), which provides the 
Bjorken scaling, comes from the probability of having the qq fluctuation of the highly virtual 
photon. There is a finite, and also scaling, contribution to ctj' from the region of < l/Q^: 

Aarip' < l/Q^ X, Q2) oc / ^a(p) oc / dp' as{p)L{p) oc asiQ')L{Q') . (18) 

Jo p Jo 

This is the ~ as{Q'')/L{Q') correction to the LLA cross section ([l7|) . Notice the strong ordering 
in the LLA cross section: 

1/Q' < p' <l/m}. (19) 

The QCD scaling violations are (logarithmically) dominated by p' ~ l/<5^- Similar analysis 

gives the longitudinal cross section 

ri r 1 /"i/™? dn^ 

aUl*N,x,Q')= dz d'p\^L{z,p)\'a{p)^— / ^^(p) 
Jo J Ji/Q^ p 

^TtJ , , -^as{p)\og[l/as{p)] oc —asmL{Q'), (20) 
which is completely dominated by p^ ~ l/Q^ (for the more discussions on this point see [36]). 

2.4 Diffraction excitation of the qq Fock state of the photon and the 'valence' qq 
component of the pomeron 

The shape of the mass spectrum from the diffraction excitation of the qq Fock state of the 
photon (Fig.3) can be estimated undoing the 2;-integration in Eq. (|^). Firstly, we note that the 
diffraction dissociation cross section is dominated by large p^ ~ i?^,mj^ [7]: 



dt 

. 1 2 



^ dz I d'p\mT{z,p)\'a{pf 



1 r aip) 1 /■i/'"/ 1 



Secondly, the invariant mass squared of the qq system equals 



M' = ( , (22) 



zl-z) 



where kg is the transverse momentum of the (anti) quark of the pair. For the crude estimation 
of the mass spectrum at > one can undo the z integration in ([7|,pTD making use of 
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k'^ ~ ~ nip so that ~ l/zp'^ and 



1 dM'^ , , 



which gives (for the more detailed derivation see [7]) 
daD,T{^* X) 



dW^dt 



\ j-ljm 

OC 



Notice the strong flavor dependence of the diffraction dissociation cross section, Eqs. (pT]), (p4D. 
For the excitation of heavy flavors the diffraction dissociation cross section is perturbatively 
calculable, for the excitation of light flavors it is evidently infrared- regularization dependent. 
(We shall encounter such an infrared-regularization sensitivity of the diffraction dissociation 
cross section over and over again. We shall demonstarte that, however, this infrared sensitivity 
can be reabsorbed into the normalization of the input structure function of the pomeron, which 
by itself will be proven to satisfy the GLDAP evolution.) The corresponding contribution to the 
structure function of the pomeron has the form [7] reminiscent of the valence structure function 
of the proton: 



iTaem(Ttot{pN} dtdM^ 



(xx{l-x) (25) 

i=0 

(The convention ([^) for the photon-pomeron cross section differs from that used in Ref. [7] by 
the factor (M^ + Q^)/M^.) Thefore, the diffraction excitation of the qq Fock state of the photon 
can be associated with DIS on the 'valence' qq component of the pomeron (see also Ref. [5]). 

2.5 Rising structure functions and higher Fock states of the photon 

The diagrams of Fig. 2 can also be reinterpreted as the Bethe-Heitler production of the qq pair 
by the photon-gluon fusion '~f*g qq. The conventional Weizsacker- Williams formula for this 
Bethe-Heitler cross section reads 

Nf ^ 

Aatoth'N^W^Q') =J2 [ dygiy,Q')a{^*g qfqf,yW') , (26) 

where g{y,Q'^) is the distribution function of the physical, transverse, gluons in the proton. By 
the kinematics of DIS 

_ {kq) _ (M^ + Q-) 
^ (pq) {W^ + Q2) • ^ ^ 
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The diagram of Fig. 2 describes the driving term of the gluon distribution in the proton, which 
here is assumed to be entirely of the radiative origin ([35,37], see also [39]). At a; ^ 1 we have 
[35,37] 

Here it is worthwhile to emphasize, that the flux of soft gluons depends only on the color charge, 
but neither the spin nor the helicity, of the source of gluons. The vertex function V{k) in the 



Aatot{j*N, W\ Q') oc J^^^ ^^-^^(7*^? - Q^, M') . (29) 



integrand of (|2q) is precisely the same as in eq. ([TT|) and describes the destructive interference 
of the radiation of gluons by different quarks bound in the color-singlet nucleon. In terms of the 
mass M of the excited qq pair, Eqs. (^),(^) give 

Because of the spin-| exchange in the t-channel, the cross section of the photon-gluon fusion 
subprocess decreases at large M^, 

^(7-«^»MV^log(«l±|^^), (30) 

the integral ( p9D converges at finite and gives the constant photoabsorption cross section at 
small X. For the closely related reason, one finds the rapidly convergent mass spectrum (^4]). 
On the other hand, if the qqg final state is produced in the photon-gluon fusion, then because of 
the spin-1 gluon exchange in the t-channel (j{'^*g —* qqg, M^) oc const, which leads to the rising 
oc log(l/x) contribution to the photoabsorption cross section [41]. Notice, that excitation of the 
qqg final state in the photon-gluon fusion can be reinterepreted as the scattering on the nucleon 
of the qqg Fock state of the photon. Therefore, one has to study the effect of higher Fock states 
of the photon. 



3 Interactions of the qqg Fock state of the photon 

3.1 Interaction cross section for the qqg state 

One can easily write down the interaction cross section for the color-singlet three-parton qqg 
state using only the color gauge invariance considerations (the separation of partons in the 
impact-parameter space is shown in Fig. 4): 

asir,R,p) = ^[a{R) + aip)] - ^a(r) . (31) 
13 



Indeed, when separation of the quark and of the antiquark is small, r p ^ R, then the qq pair 
will be indistinguishable from the pointlike colour-octet charge. In this limit 

as{0,p,p) = ^a{p) . (32) 

where 9/4 is the familiar ratio of the octet and triplet strong couplings. In the opposite lim- 
iting cases of i? = or p = the gluon and the (anti) quark with vanishing separation are 
indistiguishable from the pointlike (anti) quark and 

0-3 (r, 0, r) = a^ir, r, 0) = a{r) . (33) 

The formal derivation goes as follows: In the integrand of the cross section ([lID, the two prop- 
agators l/{k'^ + Pg) correspond to the Fourier transforms U{k) = J d^fex]i{—ikr) exp{—pGf^)/r 
and U{—k) of the (infrared-regulated) gluonic Coulomb potential. If the color charge is located 
at the position c, then U{k^ c) = U{k) exp{—ikc). Whenever the two exchanged gluons couple to 
the same parton, one gets the square of the corresponding strong charge. If the gluons couple to 
the two partons located at points fi and r2, the corresponding contribution acquires the extra 
phase factor exp[ik{f2 — fi)]. Precisely this is the origin of the factor [1 — exp{ikf)] in the inte- 
grand of eq. (|ri|) . The accurate calculation of the colour traces for the different couplings of the 
two exchanged gluons to the quark, antiquark and gluon of the qqg Fock state leads precisely to 
the cross section (|3T|). 



If $i(r, R, p, z, Zg) is the wave function of the qqg Fock state, then the corresponding contri- 
bution to <Jtot{l*p) shall equal 

^crtot{qqg,x,Q'^) = J dz cffdzg cfp \^i{r, R, p, z, Zg)\'^a{r, R, p) . (34) 

The gluon of the qqg Fock state is generated radiatively from the primary qq Fock state (Fig.5), 
and this radiation simultaneously renormalizes the weight of the qq component of the photon. 
If ng{z,f) is the number of gluons in the qqg state with the q-q separation f (we suppress the 
subscripts T and L in the |\I^(z,r)p) , defined by 

J dzgd^p\<^i{r,R,p,z,Zg)\'^ = ng{z,f)\<i/{z,r)W. (35) 

then the wave function of the radiationless qq component of the photon shall renormalize as 

\^gg{z,r)\' = \^{z,r)\'[l-ng{z,f)] . (36) 
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It is convenient to introduce 



Aa{r, R, p) = a^^r, R, p) - a{r) 



9 
8 



HR) + a{p) - a{r)] , 



(37) 



which shows how much the interaction cross section of the qqg Fock state is different from the 
cross section for the qq state. Then, atot{l*N) from interactions of the qq and qqg Fock states 
of the photon takes the form 



Up to now we have manipulated the formally divergent quantity ng{z,r) as if it were finite. 
As a matter of fact, the renormalization (^) of the radiationless qq Fock state corresponds to 
the introduction of the so-called regularized splitting functions [13] in the GLDAP evolution 
equations, and takes care of the virtual radiative corrections (a very detailed discussion of the 
interplay of the virtual and real radiative corrections and of the emergence of the running strong 
coupling was given by Dokshitzer [12], see also the review paper [40], and needs not be repeated 
here). Of course, the final result for the physical cross section, the last line of eq. (|38D, does not 
contain any divergences. 

3.2 Wave function of the qqg state and the rising photoabsorption cross section 

We are interested in the oc log(l/a;) component of the increase of the photoabsorption cross 
section 



This log(l/a;) comes from the dzg/zg integration in the domain x < Zg < 1 and we must 
concentrate on Zg <^ 1. One should not confuse the Sudakov variable Zg which is the fraction 
of the light-cone momentum of the photon carried by the gluon, with y, which is the fraction 
of the light-cone momentum of the proton carried by the same gluon. The two quantities are 
related by 




(38) 



Aailli^*N,x,Q^) = I dzd^fdzg£p\^,{r,R,p,z,Zg)\^Aa{r,R,p) . 



(39) 



ZgH = X. 



(40) 



15 



Only the diagrams of Fig. 6a in which the exchanged gluons couple to the gluon of the qqg 
Fock state give rise to Aatot{'y*N,x,Q'^) oc log(l/a;). The corresponding wave function can be 
reconstructed from the number of gluons Ug in the qqg state, which on the one hand equals 



Ug 



J dZgcfpdzcff\^i{f,R,p,z,Zg)f (41) 



and, on the other hand, can be evaluated from the Weizsacker- Williams formula 



n„ = - . - dzd'f |^(f, z)Y- / ^ ' ' [1 - exp(^A:„f)] . (42) 

Here we have used the fact that for the qq source of gluons 

Vfy = J dzd^r |^(r,2)ni - e^^iikgV)] , (43) 

and the factor 2/3 accounts for the 2 (anti) quarks in the qq state compared to the 3 quarks 
in the proton. Transformation of Eq. ( ^2]) into the configuration-space integral can easily be 
performed making use of [24] 

d^kP 



and 



which yields 



,r, ^ exp(ifcf) = / d^pplKi{pGp)Ki{pGR)^ . (45) 
{k^ + PqY J -Kp 



n. 



' 3n^J Zg 



"^'^ d'p [ dzd'f\^if,z)\'asir)pl\K,ipGP)- - K.ipcR)^^ , (46) 
J p ri 



SO that 



\^,{f,RJ,z,Zg)\' = j^^^\^{z,r)\^as{r)pl\^^^^ (47) 

The wave function (|47|) has the 1/ Zg behavior needed for the oc log(l/x) rise of the cross section. 
The color gauge invariance property of the wave function ( ^7|) is noteworthy: because of cancel- 
lations of the color charges of the quark and antiquark in the color singlet state it vanishes when 
{R — p) 0. It counts only physical, transverse gluons. For those reasons and because of the 
related color gauge invariance properties of a{r,R,p), introduction of the infared regularization 
and the modelling of the confinement by the effective mass of gluons exchanged in the t-channel 
does not conflict the color gauge invariance. 
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In the DIS the leading contribution to Aa[ll{'y*N,x,Q'^) comes from the LLA ordering of 



sizes 

1 



-r^ < < p2 ^ i?^, ^ , (48) 

V Pg 



when (here the angular averaging is understood) 

fil\K,ipGp)-^ - K^i^^GR)^\' = y , (49) 
which produces the factorized wave function 



1 .,2 4 



\^,{r,R,p)\' = -\^{z,r)\'—as{r)-. (50) 

Zg 6TT p 



Naturally, the LLA wave function (|50D does not depend on the infrared regularization parameter 
Pg- In the LLA 

Aa(r,/?,p) = S(p) = ^a(p) (51) 
and the increase of the total cross section can be written as 

AaUl*N,x,Q')= fdzd'f\^{z,r)\'a,{ry^ / ^S(p) 

J 3tT Jx Zg Jr^ p* 

= / dzd'fMz,r)\'as{r)C^r'^-^\og p' a s{p) L{p) 

= |cizc/Vlv^(z,r)|V(r)^-iL(r)-ilogg)oci^L(Q)=^logQ . (52) 

Here we have made an explicit use of the small-r behaviour of o"(r), Eq. (p^. This is the 
first instance when we encounter the expansion parameter of the Double-Leading-Logarithm 
Approximation (DLLA) [12,40,42] 

e(x,r) = ^L(r)logf-') . (53) 



We have one power of L{Q) per QCD loop (which a posteriori justifies LLA ordering (^8])) and 
one power of log(l/x) per gluon in the Fock state of the photon. 

3.3 The triple-pomeron asymptotics of the mass spectrum of diffraction dissocia- 
tion 

Our starting point is the generic formula (|^). Repeating the considerations of Section 3.2, we 
can write 

daoii* qq + qqg) 

t=0 " ^^^^.^ J^Q 



167r 



dt 



16vr fdM^ d^r -^ qq + qqg) 
J dtdM^ 
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dzd'^f r)p[l - ng{f)]a{ry + J dzdPrdZgdPp |<l>i(r, R, p, z, Zg)\^[a{r) + Acr(r, i?, p)]^ = 
(i2rfV|^(z,r)|V(r)2+ / dzd^rdzgd^p\^i{r,R,p,z,Zg)\'^[^a{r,R,pf + 2a{r)^a{r,R,p)]. (54) 



The first term in the last line of Eq. (|5^) describes the diffraction excitation of the qq Fock 
states into the low masses ~ Q^, see Eq. (^). The second term gives rise to the mass 
spectrum, which can be seen as follows: The invariant mass squared of the qqg state equals 



mj + kg ^ m?f + kl 



+ 



(55) 



Anticipating the final results, we note that the leading contribution to the diffraction dissociation 
cross section comes from the slightly modified LLA ordering 



R2 



(56) 



i.e., from 3> ^ ~ Pq. Therefore, the excitation of masses ^ only comes 

from Zg <^ z < 1, and the dzg integration in (|54D can easily be transformed into the dM'^ 
integration (see Eqs. (p7|j40|) ): 

dM'^ _dy _dzg_ 
M^ + Q^~ y ~ Zg ' ^ ^ 

where now y is the fraction of proton's momentum carried by the pomeron. The wave func- 
tion ( ^71) has precisely the needed oc 1/zg behaviour and leads to (in view of (|SB|) the term 
oc cr(r)S(r, _R, p) in ([5^ ) can be neglected) 



dan 



dtdM^ 



t=o 



dzd^r \^!{z,r)\^as{r)r'^CN 



.J_.J_.±. 

Cn 167r StT Jr2 



n 2 



F{PGP) 



(58) 



(59) 



where F[^gP) is defined by the shght generahzation of ( ^91) to allow for ^cp ~ 1: 

Pg\Ki{pgp)- - Ki{pgR)-b? = —^iPGp) ■ 
P R p* 

The form factor F{x) satisfies F{0) = 1 and F{x) oc exp(— 2x) at a; > 1. 

Firstly, we notice that the diffraction dissociation cross section depends on the infrared 

regularization, since the p integration in (|58D is essentially fiat and is dominated by large p ~ 

Rn,1/Pg- Then can take p^ = for the lower limit of intergation, and dan/dtdM"^ Eq. ( |55| ) 

factorizes into the Q^-independent dimensional constant 



A* 
^3IP 



Cn 127r2 



9x2 

4 



dp^ 



<P) 



n 2 



p- 



F{PGP) 



(60) 
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and the cross section 

a*(Q2) = I dz(fr\^iz,r)\^CNr^as{r) , (61) 

which is nearly identical to aT{'y*N, x, Q"^) Eq. ([T7|), being short of L(r) in the integrand. There- 
fore, this driving term of the triple-pomeron component of the diffraction dissociation of virtual 
photons satisfies the approximate factorization reminiscent of factorization properties of the 
triple-pomeron diagram of the conventional Regge theory Fig. 7a, 



aTii*N,x,Q'^) dtdM^ 
To the considered lowest order in the perturbation theory, the quantity A^jpi^Q'^) does not depend 
on X. A^jp^Q^) is the dimensionfull quantity, and as such it could have had a strong Q^- 
dependence, A-^jp^Q"^) ~ l/Q^ being a plausible guess if l/y/Q^ were the only scale relevant to 
deep inelastic scatering. The fact that cr*(Q^) and (7tot{l*N,x,Q^) are nearly identical, proves 
that this is not the case. Furthermore, the r.h.s of Eq. (16^) has a very smooth extrapolation down 
to the real photoproduction limit = 0, confirming the earher suggestions [26,7] that A^jp^Q"^) 
is close to A^jp{0) as measured in the real photoproduction (the more detailed comparison of 
AsipiQ"^) with the A^jp for the diffraction dissociation of protons, pions and the real photons will 
be presented elsewhere [43]). For the order of magnitude estimation of A^jp, in the dominant 
region of integration in (|60D we can take J^{^gP) ~ exp(— 2/iG'P), L{p) ~ 1 and a{p)/p^ ~ 
Cn(^s{p = with the result 

27Cn asi^fil? asi^pl? ..o^ 
^^"^ ~ 256vr2 pl = -Ji^ ■ ^^^^ 

With fiG ~ OAGeV this gives A^jp ~ 0.3 {GeV)~~^ . The experimental data on the diffraction 
dissociation of the real photons give v43ip(0) ~ O.IQ {GeV)~^ [27]. This dimensional soupling 
Agjp ^ A3ip(0) emerges as the principal nortmalization factor of the diffraction dissociation cross 
section, and Eq. (^) is a starting point of the derivation of the factorization representation (|89|) 
to all orders in the perturbation theory. 

Combining Eq. (p2| ) with the definition Eq. ( pS] ) we find the corresponding contribution to 
the structure function of the pomeron at x = Q^/(Q^ + M^) <^ 1 of the form [7] 
pPM') = , cr*{x,QMl^ ^ 16vrA3,.(0) ^(^) ^ O.OSFf ^(x, Q^) . (64) 



It describes DIS on the qq 'sea' of the pomeron. The relationship ( p^ shows a deep connection 
between the triple-pomeron component of the mass spectrum and the sea structure function of 
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the pomeron. Notice the difference between the diffractive excitation of the qq state, Fig. 3, and 
of the qqg state, Fig. 7b : in the former the pomeron couples to (anti)quarks and the DIS probes 
the 'valence' qq structure of the pomeron, in the latter the pomeron couples to gluons, and the 
DIS probes the 'sea' of the pomeron, which can be treated as having been generated from the 
'valence' (constituent) gluons of the pomeron. The diffraction dissociation cross section ( |58D and 
the pomeron structure function ( |64D are sensitive to the infrared regularization, and the normal- 
ization of both quantities contains the new dimensional parameter At^jp. The important result of 
the above analysis is that this dimensional parameter A*^jp can approximately be inferred from 
the real photoproduction data. Notice a close similarity between Eq. (|60|) for the normalization 
of the triple-pomeron mass spectrum and Eq. (^) for the normalization of the mass spectrum 
for the excitation of the qq state. However, whereas Eqs. ( ]60| , |6^ ) predict the flavor-independent 
relation between the proton and pomeron structure functions, the valence qq structure function 
of the pomeron has a strong flavor dependence [7]. Now we shall study how these conclusions 
shall change when higher order effects and QCD evolution are included. 

4 Rising structure functions and higher order Fock states of the pho- 
ton 

Generalization of an analysis of Section 3 to interactions of the higher qqgi...gn Fock states of 
the photon is straightforward. The strong DLLA ordering of gluons 

X < < Zn-i < ... < 2:1 < 2; < 1 , (65) 

^<r2<p2^....<p2 (66) 

is required to have maximal possible powers of log(l/x) and L{Q) [42,12]. To the DLLA the 
quark-loop insertions in the generalized ladder diagrams of Fig. 8 can be neglected. By virtue of 
the ordering of sizes ( p5| ) the qqgi...gn Fock state interacts like the colour- singlet octet-octet state, 
with the inner subsystem qqgi...gn-i acting like the pointlike colour-octet charge. Henceforth, 
the generalization of Eq. (|5lD is 

9 

Ao-(r, pi, p„) = S(p„) = -o-(p„) . (67) 
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The DLLA wave function is a straightforward generahzation of the wave function for 
the qqg Fock state: 



|^(r,pi,...,p„)p 



2 1 . ^ 4 r2 1 3 p2 13, .p^ 

TT^^ 2;„7r^ p4 



2;i Stt^ pf Z2 



(68) 



Notice, that the first gluon is radiated by the triplet-antitriplet colour-singlet state. The subse- 
quent gluons are radiated by the octet-octet colour-singlet state, which brings in the ratio 9/4 
of the octet and triplet strong couplings. The corresponding increase of the total cross section 
equals (here we make an explicit use of eq. (|T^)) 

4 /3^"-i 



AaS(7*iV,x,Q2 
dpi 



3tt Vtt 



J dz d'^f\'i!(z,r)\'^asir)r'^ 



-ri / \ I dpi . . f dpl^. , dzi /"^i dz2 /"^"-i dzn 
— «s(Pi) l—as{p2)-.. 



r2 Pi 



npn) ■ — — ... — 



dpi 

r2 PI 



f dfP' 

2 «5(Pl) ^as{p2)-- 
•'Pi P2 



pLi Pn ' ' "'^ ^1 "'^ ^2 

dz rfV|^(z, r)pa5(r)CArr^ 



P„-i Pn 



^asipn)L{p. 



X Zi J X Z2 



r'|\E'('z, r)| air 



^{x/z,rY 



[n + l)!n! 

Therefore, the total photoabsorption cross section can be represented as 



(69) 



dzd'^f\<if{z,r)\^a{-,r) . 



where 



a[x,r) = cr r 



E 



^(x,r) 



(70) 



(71) 



to (n + l)!n! 

is the energy-dependent dipole cross section, which generalizes the Low-Nussinov pomeron to 
the DLLA pomeron. 

The sum in (^Tp can be evaluated as (we neglect the slowly varying pre-exponential factor) 



E 



^ E 7^ « ^M^O 



leading to 



a[x, r) oc a[r) exp 



48 



as(r 



log - 

x 



(72) 



(73) 
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and 



")2 



F^iGLDAP, X, Q2) = .^_atot{x, Q^) oc exp 



48 

^0 



as(Q'' 



log (i) . (74) 



The representation (^) for the photoabsorption cross section in terms of the DLLA dipole 
cross section (ITT]) is a new result and is presented here for the first time. However, since the 



perturbative expansion (|69|) is completely equivalent to the GLDAP evolution equation for the 
structure function [11,12], the result (|74D for the DLLA growth of the structure function is 
identical to the one derived from the GLDAP evolution equations [12,40], where it appeard as a 
rising density g{GLDAP, x, Q"^) of gluons in the proton. In our light-cone s-channel approach it 
comes from interactions of the higher qqgi...gn Fock states of the photon and can be described 
in terms of the rising DLLA pomeron cross section ( |7TD for the color dipole. As a matter of fact, 
it can easily be shown that in the DLLA [30] 

a{x, r) = —r'^as{r)xg{x, = — ) (75) 

o T 

We note in passing, that the cross section (0,^,^) obviously does not satisfy the factorization 
relations usually assumed for the pomeron in the standard Regge phenomenology. The QCD 
evolution analysis of deep inelastic scattering needs not to assume any factorization of F2{x, Q^). 

5 The structure function of the pomeron and constituent gluons of 
the pomeron 

Now we consider diffraction excitation of the higher order qqgi----gn+2 Fock states of the photon 
(Fig. 7c). Large masses M of the excited state 



m 



j + kl m} + k1 k? 



^2 ^ -1 + — (76) 



Zq Zq -^-^ Zi 



will be dominated by the softest gluon contribution. The fraction y of proton's momentum 
carried away by the pomeron is related to Zn+2 as (cf. Eq. (^OD) 



Zn+2 = - (77) 

y 



and 

dM"^ dy dzn+2 



(78) 
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Using the wave function (p^), and repeating the considerations which have lead to eq. (0), we 
find the contribution of excitation of the qqgi---gn+2 Fock state to the mass spectrum of the 
diffraction dissociation and to the photon-pomeron cross section 



167r 



1 / 4 



3,n+l 



daoil* qqgi---gn+2) 



dpi , 



r2 Pi 



(^totipp) V37r/ Vtt 



dz d^f\^!{z^ r)\^as{r)r'^ 



pi Pn+l 



(^totiPP) Vvr 



dz d^r r)pa;5(r)r^ 



'"+1^ / . ^^S(Pn+2)^-^(PGPn+2) 







I 


'x Zi 


-/x Z2 





pI+1 Pn+2 
" dZn+1 



1 (i^;,, 
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a5(Pn 



+1; 



X Z\ 

dZn+1 



'1 (i2;„ 



dpi I N 
— ttslPi)-- 



r2 Pi 



dpi 

Pi-. Pn 



'n+2 



a{pn+2f^{PGPn+2) ■ (79) 



Comparison with Eqs. ([TT| ) , (|69| ) shows that the last line of Eq. (|79D can be reinterpreted as the 
dipole cross section for interaction with the pomeron treated as the two-gluon state with the 
wave function 



81 1 



1 



a(r) 



n 2 



where 



Stt'' Zg atotipp) 



Zn+2 



(80) 



(81) 



is the fraction of pomeron's momentum carried by the gluon. 

Indeed, making use of Eq. (|^) for the vertex function of the two-body system, the dipole 
cross section o'2g{p) for the scattering on the gluon-gluon state can be written as 

dzgj d^r\^2g{z,r)\^ ^ ^^^^^ as{k')[l - expitkr)] (82) 

The factor 3/2 difference from Eqs. (|TT],|I1|) is due to the ratio 9/4 of the gluon (octet) and the 
quark (triplet) strong couplings and the ratio 2/3 of the number of the constituents gluons in the 
pomeron and the number of the constituent quarks in the proton. The series of transformations 
of the integrand of (0), 

„l/p2 ^^2 



a2g{p) = 2tx p^ a s{p) j dzg j d^r\^ 

= 2tx p^ as{p) I dzg I d^f\^2g{zg,r)\^ j ^asip 
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2'np^as{p) I ^asipi) [ dzg [ d^f\^2gizg,r)\^ , (83) 
Jp2 Pi J Jp'i 



and comparison with the last hne of Eq. ( [79D complete the derivation of the wave function (^Ol). 
This is one of the central results of the present paper. 

The wave function '^2g{,Zg,r) gives the distribution of the 'valence' gluons gw{,Zg) in the 
pomeron 

gjp{zG) = J d'r\^2g{zg,r)\\ (84) 

The perturbative expansion (|79D describes the QCD evolution of this 'valence' gluon distribution. 
The above derivation holds at 2;^ <^ 1 , the region of ~ 1 requires special consideration. Only 
<^ 1 is important in the DLLA. The radius of the pomeron Rjp ~ Rn, l/po and is controlled 
by both the form factor J-'{fiGi^) and the behavior of a{r) in the saturation regime. The absolute 
normalization of the flux of soft gluons in the pomeron is given by the familiar coupling A^jp: 



Zggjp{zg) = ^ ] — T / dr' 



a(r) 



HPor) = ^ ■ ~ 0.06 . (85) 

9 (Ttotypp) 



Extrapolation of (^) also to large Zg gives an estimate of the gluon momentum integral for the 
pomeron 

IP = / dxxgjp[x) ^ — ^ ~ 0.06 . 86 

The momentum integral for the 'valence' (anti)quarks of the pomeron, discussed in Section 2.4, 
was estimated in Ref. [7] with the result {xv{qq))-sp ~ 0.1. Eq. (0) gives a few per cent estimate 
for the momentum integral for the sea (anti)quarks. Here we just emphasize that the pomeron 
needs not be regarded as a particle and on the purely theoretical grounds there are no reasons 
why the momentum integrals for gluons and (anti) quarks in the pomeron must add to 100% 
[5,7]. 

Henceforth, we have identified the three components of the input for the QCD evolution of 
the pomeron structure function: i) the valence quark-antiquark component with the structure 
function (^) (for the detailed analysis see Ref. [7]), ii) the valence gluon distribution with the 
structure function (BH), iii) the sea (anti)quark distribution given by Eq. (p3). All these input 



parton distributions are sensitive to the infrared regularization. There is nothing wrong with 
this sensitivity: the infrared sensitivity of the parton distributions is inherent to the QCD- 
improved parton model. In the conventional parton model phenomenology it is hidden in the 
parametrization of the parton densities at small factorization scale, which then is used as an 
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input in the QCD evolution analysis of the scaling violations. The important finding is that the 
absolute normalization of the sea and gluon distributions in the pomeron is determined by one 
and the same fiavor-independent dimensional constant Agjp which must approximately be equal 
to the triple-pomeron coupling as measured in the real photoproduction. The normalization 
of the valence qq structure function of the pomeron is given by the very similar but fiavor- 



dependent dimensoinal constant (cf. equations ( pi]) and (|58| , |60D ) . In the above DLL A analysis 
we omitted the quark loops in the ladder diagrams for the pomerons. These quark loops shall 
automatically be included in the GLDAP-evolution calculation of the structure function of the 
pomeron starting with the above described input parton distributions in the pomeron. 

6 Flux of the QCD pomerons in the proton 

To complete our analysis we must replace the Low-Nussinov two-gluon pomeron in the lower 
part of diagrams of Fig. 7c by the full QCD pomeron - the sum of the triple-ladder diagrams 
of Fig. 7d. This is done by replacing the dipole cross section cr(p) by cr{y,p) in the last line of 
Eq. (|79|), where y is a fraction of the proton's momentum carried by the pomeron. Then, the so 
calculated diffraction dissociation cross section will differ from that of Section 5 only by the y 
dependent factor 

My) = Idp'Hy^P)/p'?nf^cP) (87) 
Idp^cr{p)/p^fJ-{^^Gp) ^ ^ 

What is the proper interpretatation of /ip(y) ? 

We would like to preserve the most important result of the above analysis - the repre- 
sentation of the diffraction dissociation cross section through the GLDAP-evolving structure 
function of the pomeron. The scaling variable of the photon-pomeron scattering equals 
xjp = Q^/{Q^ + M^) = x/y, so that 

dM^ dy 



M'^ + Q^ y 

With allowance for the factor fjp{y) Eq. (87) for the diffraction dissociation cross section can be 



written down in the factorized form 

dao 



dt dy 



t=o 167r Q"^ y ^ y 
which has the conventional parton model representation with f-sp{y)/y being the flux of pomerons 
treated as partons of the proton. In order not to introduce any spurious dependence on the 
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kinematical variables x and y, the coefficient atot{pp)/^Q'^ in ( P^ must be taken constant, for 
instance fixing <Jtot{pp) = 40 mb. Because the pomeron is not the particle with the well defined 
couplings (residues) and spin, and because in QCD the pomeron does not factorize, the regge- 
theory convention (|1]) is not unique. The coefficient atot{pp) / ^Q^c is the convention-dependent 
normalization constant for the correct dimensionality of the diffraction dissociation cross section 
in terms of the dimensionless structure function or vice versa. The absolute normalizations of the 
flux of pomerons and of the pomeron structure function too are the convention dependent ones: 
it is always the product of the two quantities which enters the experimentally observable cross 
sections. Eq. (^9]) shows how the QCD evolution effects and the rising dipole cross section cr(x, p) 
modify the law @ for the mass spectrum in the triple-pomeron region. The factorization 

( |8^) bears certain semblance of the usual Regge-theory factorization in the triple- Regge region. 
We emphasize that we have derived (pOf ) neither assuming nor using any of the Regge-theory 
factorization relations. 

The three pomerons in the triple-pomeron diagram are described by different QCD ladder 
diagrams. The top pomeron of Fig. 7d is in the LLA regime: the relevant sizes vary along 
the ladder from ~ bottom cell of the ladder down to ~ in the top, 

quark-antiquark, cell of the ladder. For this reason we can introduce the structure function of 
the pomeron. By contrast, the exchanged pomerons in the lower part of Fig. 7d are the soft 
pomerons: since p^^2 ~ ~ ^^e sitiuation is reminiscent of the pomeron contribution to 
the typical hadronic cross section, see eq. (0). The predictive power of QCD for the hadronic 
total cross section is still very limited [44,19,17,45]. (We shall comment more on Lipatov's 
work on the QCD pomeron [29] below.) The empirical observation is that the hadronic cross 
sections and the real photoabsorption cross section [46] have a very weak dependence on energy 
u ~ rriN/x, much weaker compared to a steep rise of the DIS structure functions with 1/x. The 
dipole cross section ([Til) is consistent with this property: it is essentially flat vs. 1/x at large, 
hadronic, size r, and the smaller is the size r, the steeper is the rise of a{x,r) with 1/x. 

This rise has a certain impact on the radius of the pomeron. Namely, the replacement of 
a{r) by a{x,r) leads to the effective wave function of the pomeron 



|^ip(y,2G,r) 



2 



•^(/icr) . (90) 



Svr^ Zg (TtotiPP) 

With the rising dipole cross section a{y,r) Eq. (4.7) the ratio a{y,r)/r'^ will rise towards small 
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r, so that the effective radius of the pomeron Rip{y) will decrease as ?/ — 0. The GLDAP 
evolution equations hold at R^Q^ ^ 1, where Ri are the relevant hadronic radii. If one would 
like to formulate the input for the GLDAP evolution at certain fixed then, because of the 
decreasing radius of the pomeron, the GLDAP applicability condition will be violated at very 
small values of y < ydQl) such that 

RwivciQ'fQ' = I . (91) 

The factorization of the diffraction dissociation cross section into the flux of pomerons and the 
structure function of pomerons Eq. (pQ]) will break down a.t y < yc{Q'^). Our criterion (pT]) for 



the breaking of the GLDAP evolution is different from GLR criterion ([19], for the review and 
references see [22]) and its implications will be studied in the subsequent publications. 

A brief comment on the target dependence is in order. In the Regge theory, in the approxi- 
mation of exchange by the single factorizing pomeron. 



1 daoil* X +p) 



f=o '^toi(7rvr) dtdM"^ 



(92) 

t=o 



CTtotipp) dtdM"^ 

In our s-channel approach the target-dependence enters through the target- dependent dipole 
cross section ([TT]) , which is roughly proportional to the number of quarks in in the target hadron b 
(To the extent that the quark-quark separation in the proton and the antiquark-quark separation 
in the pion are similar, the Low-Nussinov cross section (^) reproduces the additive quark model 
[33,14,15].): 

(j{{qq)p,p) ^ (JtotiPP) ^ o-tot(p7r) ' 
so that despite the lack of factorization, the QCD pomeron roughly satisfies the relation (|9^ 



7 Unitarization of the rising structre functions 
7.1 Rising cross sections and the s-channel unitarity 

The rising cross section a{x, r) Eq. ( |7TD conflicts the s-channel unitarity at sufficiently large 1/x. 
The s-channel unitarity constraint is best formulated in the impact-parameter representation 
(the partial-wave expansion) and reads 

T{b) < 1 . (94) 
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The profile function T{b) is related to the elastic scattering amplitude f{q) such that 

f{q) =2ij d^bTib) exp{-iqb) = latot exv{-]^B,if) . (95) 

Here Bei is the diffraction slope of the elastic scattering. The Gaussian parametrization (p5|) is 
viable for the purposes of the present discussion [47] and gives 

The profile function of the bare pomeron excahnge ro(6) defined for the rising cross section ([7T|) 
will overshoot the s-channel unitarity bound at a sufficiently small x. 

The are no unique prescriptions how to impose the s-channel unitarity constraint on the 
rising cross sections. The often used procedures are the eikonal [48,49] 

m = 1 - exp[-ro(6)] = E ^-^^^oiby (97) 



and the /C-matrix [50,44] 



s-channel unit arizat ions. Both produce T{b) which satisfies the unitarity bound (plf). To the 
leading order in the s-channel unit arizat ion, the unitarized profile function reads 

m^ub)-lxUbr (99) 

with X = 1 foi' the eikonal unitarization, and x = 2 for the /C-matrix unitarization. The eikonal 
unitarization is being routinely used in high-energy physics [51] and sums the s-channel iterations 
of the bare pomeron exchange (Fig. 9a) when only the elastic scattering intermediate states are 
included in the s-channel. Besides the elastic scattering states as in Fig. 9a, one must include 
the inelastic intermediate states of Fig. 9b, which correspond to the diffraction dissociation of 
the target nucleon. These inelastic intermediate states lead to an enhancement of the double 
and higher order rescattering terms in expansions (p7| , |98D [52,49,53]. If one starts with the 
eikonal unitarization (which is an assumption), and includes the corrections for the diffraction 
dissociation of the target nucleons, then [53] 



In fact, the /C-matrix prescription (^) was obtained in [44] starting with the eikonal unitarization 
of ttN scattering and including the inelastic intermediate states in the QCD inspired model of 
the diffraction dissociation of pions. 
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7.2 Shadowing correction to the proton structure function 

Now we shall discuss the unitarization (shadowing, absorption) effects in DIS, taking full ad- 
vantage of the diagonalization of the S-matrix in the (p, 2;)-representation, which allows us to 
impose the s-channel unitarization on all the multiparton cross sections (jn(r, pi, p„) at all 
the values of r and pj. Although the unique s-channel unitarization procedure is lacking, we 
can still develope a sound phenomenology. We identify the cross section ( fflD which leads to 
the GLDAP-evolving structure function F2^\gLDAP, x, Q"^) with the bare-pomeron exchange. 
The bare-pomeron structure function F2^\gLDAP, x, Q'^) is the linear functional of the density 
of partons in the proton: 

F^''){GLDAP,x,Q^) = -^atotiGLDAP,x,Q^) 



e^iX \qi{GLDAP, x, Q^) + q^{GLDAP, x, 



flOll 



The construction of the unitarized photoabsorption cross section goes as follows. For each Fock 
state we define the bare ro(&) and the unitarized r{b) profile functions and the bare cxo and the 
unitarized cross section a^^^: 

a(^) = 2 J d%T{h) = 2 j (fbVoib) -2 J (fb[To{b) - T{b)] = ctq - Aa^'^^ (102) 

Let us derive the shadowing (unitarity) correction to the scattering of the qq Fock state of 
the photon. To the leading order in To{b), Eqs. (p9D,(p6D give 



and the shadowing correction to the total photoabsorption cross section equals 

Aatt\l*N,x,Q') = I dzd'p\^iz,p)\'Aa^^'\p) ^xj dzd'p\^>{z,p)\^^^^ 



dtdM^ Bd{M'^) dtdM^ 

dM^ atotipp) 



t=o 



Here B£,{M^) is the diffraction slope for the diffraction excitation of the mass M. The shadowing 
correction to the total photoabsorption cross section equals the diffraction dissociation cross 
section times the enhancement parameter x ~ (1 ~ 2). The generalization of eq. (|104|) to 
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interactions of the higher Fock states of the photon is straightforward. Making use of eq. 
we obtain the shadowing correction to the structure function of the proton 



AF<-'(x. Q^-) = xfiM r '^My)Fr\±, Q=)^ . (106) 



(The slope B^jp and the end-point i/m of the pomeron distribution will be defined below.) Ignore 
for a minute the mass-dependence of the slope Bd{A/P). Since F^lxjQ"^) satisfies the GLDAP 
evolution, the convolution representation (|105|) implies that the shadowing correction to the 



proton structure function also satisfies the GLDAP evolution equations! Experimentally, in all 
the hadronic reactions and in the diffraction dissociation of real photons the slope i?£)(M^) 
exhibits similar dependence on the excited mass M [27]: in the triple-pomeron region the slope 
is constant to a good approximation, 

Bd{M') = B,jp ^ \B,i{hN) , (106) 

whereas in the resonance excitation region 

Bd{M^) ~ B,i{hN) . (107) 

In the DIS the counterpart of excitation of resonances is the excitation of the qq Fock states 
of the photon, for which we expect the slope ( p.O?l ), whereas for the higher Fock states and 
heavier masses the slope ( |10(j| ) is more appropriate. These assumptions can be tested at HERA. 
Consequently, as compared to the pomeron structure function as measured in the diffraction dis- 
sociation, in the shadowing structure function (|105|) the 'valence' qq component of the pomeron 
enters with the suppression factor B^jp/ Bd^M"^) ^ 1/2, which does not effect the QCD evolution 
properties. We conclude, that the unitarized structure function of the deep inelastic scattering 

^r^(a;,<5') = F^''\gLDAP,x,Q^) - AFi'''\x,Q^) (108) 

satisfies the linear GLDAP evolution equation. 

7.3 Brief phenomenology of the shadowing correction to the proton structure func- 
tion 



According to Eq. ( |104D , the relative shadowing correction to the proton structure function equals 



the fraction wod oi DIS wich goes via diffraction dissociation of photons times the enhancement 
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parameter x- In the diffraction dissociation events the proton changes its longitudinal momentum 
Pl little, ^pl/pl = y < 2/m ^ 0.1, and appears in the final state separated from the hadronic 
debris of the photon by the rapidity gap 

Ar] = log (109) 

The standard definition of the diffraction dissociation corresponds to Ar] > Arj^in = 2 — 2.5. 
The maximal value of the rapidity gap is r^max = log(l/^)- The estimate of Wdd is particularly 
easy when the pomeron and proton structure functions are approximately constant. In this 
case fwiy) = 1, the rapidity gap distribution is fiat which is a signature of the triple-pomeron 
mechnaism [1], and combining equations (^) and ( |104[ ) we find [7] 



wdd{x)^— / — = — / dA7]=— log — , (110) 

which is roughly Q^-independent. Numerically, A3ip(0)/i?3ip ~ 0.03, and in Ref. [7] we gave an 
estimate wdd ~ 0.15 at x ~ 10~^ and ~ 30 (GeV/c)^. This prediction is consistent with the 
first determinations of wdd by the ZEUS collaboration at HERA [28]. 
For a somewhat more realistic evaluation of wdd let us assume that 

/ \ 

My) - [jj , (111) 

where A = aip(O) — 1 ~ 0.1, as suggested by the pomeron phenomenology of the hadronic cross 
sections [16,17]. We also assume that at small x the struicture functions rise oc {1/xY with 
the same exponent 6 for the proton and pomeron. The analysis of Ref. [30] gives S{Q'^) ~ 0.21 
at g2 = 4 (GeV/c)^ and 6{Q^) ~ 0.31 at = 15 (GeV/c)^. The experience with the QCD 
evolution analysis suggests that the ratio F2^\x, Q'^)/F2^\x, Q"^) only will weakly change with 
Q^, so that Eq. (|64D can be used for the relative normalization of the proton and pomeron 
structure functions. Then, 

^3IP(0) Z"'?™^ \ /\ A M ^3IP(0) 1 / ^ ^ "^1 



wdd{x) ^ — / rfA77exp[-7(Ar7 - Ar]min)] = 

-D3IP JAr]min -D3IP 7 



1- - 

\ymy 



(112) 



and to the extent that 7 = 6 — 2 A <^ 1, and 'y{f]max — Arjmin) ^ 1, we have still an approximately 



fiat rapidity gap distribution and again obtain the estimate ( |110| ) for wdd- Consequently, we 
predict rather large shadowing effect in the proton structure function 

AFi''')(x,g2^ 



x,Q2) 
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XWdd{x) , (113) 



which persists at all . In the kinematical range of the DIS at HERA the shadowing effect 
can be as large as ~ 30%. The more detailed phenomenology of the shadowing corrections is 
presented in Ref. [30]. 



7.4 Unitarization and shadowing correction to the parton densities 

Since F2^\x, Q"^) satisfies the GLDAP evolution, the shadowing correction to the proton struc- 
ture function can be reabsorbed into the modification of the parton densities in the protons. For 
instance, the shadowed density of gluons in the proton shall equal 

g{x^ Ql) = g{GLDAP, x, Ql) - T ^/ip(|/)^ip(-) • (114) 

Similarly, the valence and sea qq distributions in the pomeron will modify the sea quark distri- 
bution in the proton. The detailed phenomenology of the shadowing corrections to the parton 
distributions in the proton will be presented elsewhere [31]. Here we just notice that whereas 
the GLDAP defined parton distributions satisfy the momentum sum rule 

^ (xp) = dxx\g{GLDAP,x,Q'^) + Y,[qi{GLDAP,x,Q'^) + qi{GLDAP,x,Q'^ 



= 1 , 
(115) 

because of the shadowing correction this sum rule does not hold for the experimentally measured 
shadowed (unitarized) parton distributions. A crude estimate of violation of the momentum sum 
rule (|1T5D is 



AS,. = 1 - E i^p) ~ Z/m ■ 4^ ~ 0.003 (116) 

P='},<],9 



With the ~ 2% systematic normalization errors in the most accurate measurements of F2^\x, Q'^), 
presently the momentum sum rule can not be tested to better than 5% [55]. The concept of the 
fusion (recombination) of partons must be used with much caution. For instance, the shadowing 
correction to the density of gluons Eq. (|115|) is not proportional to g{x, Q'^Y as it is often stated 
in the literature ([19-22], see below Section 7.6). Indeed, the shadowing term is proportional 
to A^jp, the integrand of which is oc [cr(y, r)/r^]^ oc [xg^x.Q"^ = 1/r^)]^ and the integration 
is dominated by large hadronic values of r ~ Rjp and small virtualities of the fusing gluons 
Q|, ~ l/Rl> [23,24,35]. 

Because of the shadowing the parton distributions inferred from the GLDAP evolution 
analysis of the DIS structure functions will be different from the operator-product expan- 
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sion (OPE) defined parton distributions, which define the impulse approximation component 
F2^\GLDAP, X, Q^) in Eq. ( |101| ). To this end, an analogy with the comparison of the electron- 



nucleus and proton-nucleus scattering is instructive: The elastic eA scattering is described by 
the sum of the impulse approximation diagrams of Fig. 10a and is a linear functional of the 
nuclear charge density. The eA scattering amplitude measures the charge of the nucleus which 
equals the sum of charges of its constituents (nucleons). Choosing an appropriate external field, 
one can study the whole sequence of the nuclear matrix elements which will be sensitive to the 
momentum distribution of nucleons in the nucleus. For instance, considering the scattering of 
the nucleus on the gravitational center one can derive the momentum sum rule that the con- 
stituent nucleons carry the total momentum of the nucleus [35] . Under the strong condition that 
the scattering in external fields is described by the impulse approximation, i.e., by the exchnage 
of the single quantum of the external field, having measured the amplitudes of scattering in 
a variety of external fields one can reconstruct the momentum distribution of nucleons in the 
nucleus. One would recognize in the above the standard OPE definition of the parton densities 
(for instance, see the textbooks [56]). In the pA elastic scattering the impulse approximation 
amplitude /a(^ = AfN{(l)GA{(l), where Ga{(i) is the body form factor of the nucleus, has the 
profile function T^ib) ~ A^^^ which grossly overshoots the unitarity bound (0). Consequently, 
the pA scattering amplitude is subject to the large unitarization corrections (Fig. 10b) and is a 
nonlinear functional of the nuclear matter density. 

In this context, the GLDAP approach corresponds to the impulse approximation and Eq. ( |101D 



gives the linear relationship between the Compton scattering amplitide and the parton densi- 
ties. The shadowing term is an apparently nonlinear functional of the density of partons in the 
proton, but we have proven that this nonlinearity can be cast in the form of the renormalization 
of the parton densities, with retention of the linear GLDAP evolution properties. 

7.5 The higher order unitarity corrections and fusion of partons 

The higher order unitarity corrections, i.e., the multiple-pomeron exchanges Figs, lb. Id, 9, do 
technically give rise to the photon-miltipomeron interactions, which cast shadow on the very def- 
inition of the photon-pomeron cross section and the pomerons structure function Eqs. (|I|,0). The 
remarkable observation is that ond can still describe the diffraction dissociation cross section in 
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terms of the pomeron structure function and the factorization representation Eq. (PP|), and these 
7*(?7,IP) interactions only do shghtly modify fspiv) and the simple relationship ( |105| ) between 
the shadowing staructure function and the pomeron structure function. Let us start with the 
unitarization of the diffraction dissociation cross section. The s-channel iterations of the QCD 
pomeron exchange to the left and to the right of the unitarity cut in Fig. 11 do separately summ 
up to the unitarized dipole cross section. For the qq Fock state one must unitarize a{y,r), for 
the qqgi-..gn Fock states one must unitarize T,{y,r) = |cr(?/,r). Barring the qq state the flux of 
pomerons fjp^ (y) which enters the diffraction dissociation cross section must be calculated with 
the substitution 

aiy.rf Q\^^\y,rr (117) 
in the pomeron wave function (|90D , so that 



fTpiy)=[n] , ' r . N / oi2 ' , ^ • (118) 



.9/ J dr'^ [a{r)/r'^] J^{^g^) 
Apart from this minor change, the perturbative QCD expansion (ffOD will retain its form. 

Similarly, in the case of the shadowing correction. Fig. 12, the higher order unitarity correc- 
tions are accounted for by the substitution 

<r(y,r? . {^^' . AE<"-)(y,r) = 1^21 . . pfe..) - sTOfo,.)] (119) 

SO that xfipiy) in Eq. ( |105| ) shall be replaced by 

Kim ts/) .ilP ygj Jdr'l^{r)/p^-fT(pcp) 

Again, for the qq state one must unitarize cr{y,r). The higher order unitarity corrections make 
the two fluxes flp\y) and fjp^\y) slightly different both in the absolute normalization and in 
the ^/-dependence. 

Equations (|118D and (|120|) summ in a very compact form all the multiple-pomeron exchanges 
in the s-channel (Figs. 11, 12). The origin of this remarkable result is simple: the interaction 
cross section of the n— parton Fock state of the photon is dominated by the spatial extension 
of the softest gluon, which acts as a constituent gluon of the pomeron. This corresponds to 
the dominance of the multipomeron exchange diagrams of Fig. 13. Consequently, the unitariza- 
tion affects only the normalization of the pomeron wave function but not the QCD evolution 
properties of the pomeron structure function. 
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The shadowing structure function can be reinterpreted in terms of the fusion of partons 
from the overlapping pomerons emitted by the same nucleon, which reduces the total density 
of partons. The fusion of partons from different nucleons of the nucleus was first introduced in 
1975 [54] and remains a viable mechanism for the nuclear shadowing in deep inelastic scattering 
[23,24,35]. However, this interpretation must be taken with the grain of salt. The bare, GLDAP 
cross section ([75|) is a linear functional of the density of gluons, the unitarized cross sections 
contain the terms oc {—l)"''^^[xg{x,r)]"- which are sign-alternating. In more 
general terms, the multi-gluon exchange contribution is proportional to the many-gluon density 
matrix, the elements of which are not neccessarily positive defined. Evidently, this quantum- 
mechanical property is missed in the probabilistic approach to fusion. 

7.6 Unitarization and linear GLDAP versus nonlinear GLR evolution equations 

There was much discussion of the unitarization of rising structure functions in the framework 
of the so-called Gribov-Levin-Ryskin (GLR) nonlinear evolution equation ([19], for the recent 
review with many references see [22]). Here we briefly comment on the origin of the nonlinear 
term in the GLR equation, following the standard derivation of the evolution equations [11-13]. 
(We only consider a; ^ 1 of our interest.) One starts with evaluation of the derivative [12] 

^'dQ^^'^"^' Q'^[QVto,(7*Ar, X, Q')] . (121) 



In Section 2 we have decomposed Q'^atot{'J*N,x,Q'^) into the non-LLA component (|TH), which 



we can neglect, and the LLA component (p!7D, in which all the explicit dependence on is 
concentrated in the integration limit. This is the crucial point, since taking the derivative (|121|) 
and making use of Eq. ([75| ) one obtains one of the small-x GLDAP equations 

Q'^ixqix, Q')] cx gV(l/g) cx xas{Q')g{x, Q') , (122) 

in which the both r.h.s and l.h.s are evaluated at the same value of Q^. Notice, that this property 
is a result of the singular behaviour of the integrand in Eq. (p!7|). 

To the contrary, the integrand of the leading shadowing correction ( |104 .21) is a smooth 



function of p. Furthermore, it is dominated by the contribution from large p ~ Rjp. For 
this reason, it would be illegitimate to enforce the LLA limit of integration > 1/Q^ in the 



shadowing correction ( |104| ). If, nonetheless, one proceeds so, then differentiation in the first line 
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of Eq. ( |104D will give 



(123) 

The familiar form of the GLR nonlinear shadowing correction to the GLDAP equation for the 
density of gluons [19], 

(3=^|.A<,«"')(..Q')| oc ^ r ^Ji Wy,QY . (124) 

is different from our result (|105| , |89D . Evidently, neglecting the contribution to the shadowing 



term from < 1/Q^ and/or the oc 1/Q^ corrections to the leading form of the wave function 
in Eqs. (|21| , |104| ) can not be justified, which makes the GLR equation highly questionable. It is 



interesting to notice, that Mueller and Qiu [20] had already have expressed similar doubts on 
the validity of the GLR nonlinear equation (see also the recent preprint by Levin and Wiisthoff 
[21]). The GLR term is a part of the oc l/Q"^ corrections to the leading shadowing term given 
by our Eq. ( |105| ), see also the above discussion of the fusion of partons in Section 7.5. 

8 Lipatov's pomeron and diffractive deep inelastic scattering 

Lipatov and his collaborators have shown [29,42,57] that the QCD pomeron of the Regge limit 
Q^, = const, 1/x oo has the intercept 

121og2 , , 

aiP = l + Aip = l + ~ 1.5 . (125) 

TT 



which would have given a very rapidly rising total cross section 

crtot cx - . (126) 



In our scenario of diffractive DIS we have assumed slow growth of the dipole cross section ([7T| ) 
with 1/x at large, hadronic, values of r ~ R^. Firstly, this is consistent with, and supported by, 
the lack of rapid rise of xg{x,Q'^) at small (see Eq. ([75|) ) found in the latest QCD analysis 
of the parton distributions in the proton [55]. Secondly, this assumption is perfectly consistent 
with the Lipatov's theory of the pomeron. 

Lipatov's starting point is the lowest order QCD cross section for the scattering of the two 
color dipoles ri and r2 

32 r (Pk r ^ 1 r -* i 

(^ddin, ra) = —al J ^ ^ \l - exp(-zfcfi)J [l - exp(ii/cf2)J . (127) 
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In the frozen-a^ approximation, the profile function of the dipole-dipole scattering has the form 
(here the impact parameter b is defined relative to the centers of the projectile and target dipoles) 

Tddin, f2, b) cx Ko{fiG\b + s\f + Koificlb- s\f + Koificlb + A\f + Koificlb - A\f 
-Ko{fiG\b + s\)Ko{fiG\b + A|) - Mficlb + s\)Ko{fiG\b - A|) 
-Koificlb- s\)Ko{fiG\b + A|) + -Ko{fiG\b - s\)Ko{fiG\b - A|) , (128) 

where s = (ri + r2)/2 and A = {fi — r2)/2. In the limit of 



(129) 



one can use Ko{x) oc log(x). Then, the dependence on /ig in Eq. ( |128D disappears and the 
dipole-dipole profile function acquires the conformally invariant form [29] 

+ A||6- A| 



^dd{ri,f2,b) oc log 



log 



|6 + A||6- 


- A| 




fi| 


1^2 





(130) 



\b + s\\b-s\ 

It differs from (|128|) by terms which give vanishing contribution to the dipole-dipole cross section 



after the d'^b integration. In the same conformal-invariant limit, also the wave functions of the 
many body Fock states will become the scale invariant functions, see Eq. (^7]). Because of this 
scale invariance no specific size ordering dominates in the higher order QCD diagrams. Using the 
powerful technique of the conformal field theories, Lipatov has found that with allowance for the 
running QCD coupling the QCD pomeron corresponds to the series of poles, which accumulate 
at j = 1 in the complex angular momentum plane, with the intercept of the rightmost singularity 
given by eq. ( [125| ). Calculation of residues of these Lipatov's poles is as yet lacking. 

In the scenario of Ref. [44] the dominant, constant, component of the hadronic cross section 
comes from the j = 1 cross section (pT]), whereas the contribution of poles with Ajp > has 
a small residue. Indeed, the dipole cross section a{p) is given by the expectation value of the 
dipole-dipole cross section over the target nucleon state. 



crir] 



d'R,d^R2\^N{Ri, R2, R3)\'crdd{r, Ri± - R21.) , 



(131) 



and receives the dominant contribution from the separation of quarks in the proton in the 
nonscaling region of {Ri±—R2±Y ~ R'n, 1 / f^c- Lipatov's conformal pomeron with the intercept 
( |125| ) only is applicable in the scaling region of {Ri± — R2±Y ^ 1//^G5 then Eq. ( |131j ) gives 
a natural explanation for the small residue of the rightmost pole, corroborating the scenario [44] 
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(for the evaluation of the effective intercept of the pomeron in the hadronic scattering at finite 
/iG see Ref. [43]). 



Lipatov's analysis and the intercept (125) with the running coupling as{t) are expected to 
be directly applicable to the hard elastic pp, pp scattering at \t\ ^ m^, |t| ^ s. The dominance 
of the rightmost pole of the pomeron would have implied 

^^£^M oc ,2(«ipW-i) ^ ,1 . (132) 
dt ^ ' 

The experimental lack of such a rapid growth of the differential cross section of elastic scattering 
with energy is another strong evidence for the small residue of poles with Ajp > 0. 

The DLLA pomeron and Lipatov's pomeron summ (in a somewhat different kinematical 
regime) the gluon ladders with two gluons in the t-channel. They both share the unpleasant 
property of running into conflict with the s-channel unitarity, which shows that the initial 
DLLA and/or LLA selection of the subset of the perturbative QCD diagrams for the pomeron is 
unsatisfactory and open to a criticism. As a matter of fact, such a construction of the pomeron 
lacks self-consistency, since the sub-leading unitarization corrections do become larger than the 
'leading' single-pomeron exchange. From the numerical point of view, the above presented 
phenomenology seems to be viable in the kinematical range of HERA: according to our analysis 
in Section 7.3, the unitarization effects are numerically large, but still significantly smaller than 
the leading GLDAP component of the structure function (for the more detailed phenomenology 
of the shadowing see Ref. [30]). More theoretical work on unitarization is badly needed. Some 
16 years ago, Matinyan and Sedrakyan [58] have ponted out the potential significance of the 
mutiparticle Regge singularities (Fig. 14). In the recent very interesting paper [59], Bartels 
has shown that such a correlated four-gluon exchange leads to the new singularities in the 
complex angular momentum plane which are missed when the eikonal-like exchange by the two 
rightmost poles of the Lipatov's pomeron is considered. The phenomenological implications of 
this result for the shadowing corrections to the proton structure functions are not yet clear, 
since the quantitative understanding of residues of Lipatov's poles in the DIS regime is lacking. 
Our minimal-regularization approach, in which the radiative generation of the glue in protons 
([35,37], see also [39]) is closely related to the slow rise of a{x,r) at large, hadronic, values of 
r, leads to a good parameter-free description of the proton structure function statring with the 
SLAC-NMC range of x and [23,24,35] and down to the kinematical range of HERA, where 
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we have obtained [30] a good agreement with the HERA results 

onFf^(a;,g2) [60]. 

9 Conclusions and discussion of results 

Our principle conclusion is that the diffraction dissociation of virtual photons in DIS can be 
described as the DIS on pomerons with the well-defined and GLDAP-evolving structure func- 
tion. Our analysis completes the proof of the parton model phenomenology of the pomeron put 
forward some 8 years ago by Ingelman and Schlein [2]. Furthermore, we have shown that such 
a description persists beyond the single-pomeron exchange approximation. Our new result is 
that we have identified the valence qq, the valence glue and the sea qq parton distributions in 
the pomeron, which are to be used as an input in the QCD evolution of the pomeron structure 
function. We have found that the normalization of the valence glue and sea in the pomeron is 
fixed by the single dimensional coupling A^jp, which is sensitive to the infrared regularization. 
Our principle finding is that this coupling A^jp (and the corresponding triple-pomeron coupling 
AzwiQ"^) which we have shown only weakly depends on Q"^) must be approximately equal to the 
triple-pomeron coupling A3ip(0) as measured in the diffraction dissociation of the real photons 
[27]. This approximate equality Agjp fa ^3]p(0) was conjectured long ago [26] and has been a 
basis of the successful phenomenology of the nuclear shadowing in DIS [23,24,35]. This equahty 
also was used in the prediction [7] of the rate of the diffraction dissociation in DIS, which is in 
good agreement with the first data by the ZEUS collaboration [46]. Important implication of 
separation of the infrared-sensitive input structure function of the pomeron from the hard QCD 
evolution effects is that the jet activity in the DIS on the pomeron must be similar to that in 
the DIS on the proton. 

We have derived the unitarity (shadowing) correction to the proton structure function at 
small X and have demonstrated, that the unitarized structure function satisfies the conventional, 
linear, GLDAP evolution equations. We emphasize the intrinsic simplicity of our light-cone s- 
channel formalism used in this derivation. Firstly, our formalism implements in a very simple 
way the color gauge invariance constraints. Secondly, exact factorization of the photoabsorption 
cross section into the wave function and into the (multiparticle) dipole cross section allows an 
easy identification of the partial waves of the dipole cross section as an object of the s-channel 
unitarization. Thirdly, we took full advantage of the diagonalization of the scattering matrix 



39 



as a function of the transverse separation and longitudinal momenta of partons in the multi- 
parton Fock states of the photon. This enabled us to easily impose the s-channel unitarization 
on the total cross sections of all the multiparticle Fock states of the photon. This also enabled 
us to identify the constituent gluon wave function of the pomeron, which gives a very economic 
description of the shadowing process in terms of the single parameter A^jp, which is under the 
good control as it is related to the triple-pomeron coupling A^jplO) known from the real pho- 
toproduction experiments. We have shown how the multipomeron exchanges in the shadowing 
structure function and in the diffraction dissociation can be summed in a very compact form 
which only renormalizes the effective flux of pomerons in the proton. 

Applications of the above formalism to the hadronic scattering problem will be presented 
elsewhere [43]. 
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Figure captions 

Fig.l - The single- and multiple-pomeron exchange contributions to the (a,b) elastic scattering 
and (c,d) the diffraction dissociation amplitudes. 

Fig. 2 - The lowest order QCD diagrams for interaction of the qq Fock state of the photon with 
the target nucleon. In all the figures the wavy, sohd and dashed lines are for the photon, 
(anti) quarks and gluons, respectively. 

Fig. 3 - One of the 16 lowest order QCD diagrams for the inclusive cross section of the forward 

diffraction dissiociation of the qq Fock state of the photon. The vertical dashed lines shows 
the unitarity cut corresponding to the diffractively excited state. 

Fig.4 - The spatial structure of the qqg Fock state in the impact-parameter plane. 

Fig. 5 - Scattering of the qqg Fock state of the photon on the nucleon by interaction of its radia- 
tively generated gluon. 

Fig. 6 - Different couplings of the exchanged gluons to the color-octer qq pair and the gluon of 
the qqg Fock state of the photon. 

Fig. 7 The triple-pomeron diagrams for the diffraction dissociation of virtual photons in the deep 
inelastic scattering: 

[a] - The triple-pomeron diagram which describes the oc component of the mass 
spectrum in the triple-Regge phenomenology of diffraction dissociation. 

[b] - The driving term of the triple-pomeron mass spectrum in QCD - the diffraction 
excitation of the qqg Fock state of the photon. 

[c] - Diffraction excitation of the many-particle Fock states in the Low-Nussinov ap- 
proximation for the exchanged pomerons. 

[d] - The same as (c) with the exchange by the full QCD pomerons. 

Fig.8 - Gluon-ladder representation of the DLLA pomeron in QCD. 

Fig. 9 - s channel iteration of the pomeron exchange: 

[a] - In the approximation of elastic intermediate states. 
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[b] - Contribution of the inelastic intermediate states (diffraction dissociation of the 
target) to the s-channel iteration of the pomeron exchange. 

Fig. 10 - The archetype operator product expansion: 

[a] - The impulse approximation diagram for the electron-nucleus scattering which is a 
linear probe of the nuclear charge distribution. 

[b] - Multiple scattering diagrams which unitarize the proton-nucleus elastic scattering 
amplitude. 

Fig. 11 - Absorption (unitarization) corrections to the diffraction excitation of the qqgi----gn Fock 
state of the photon. The vertical dashed line shows the unitarity cut. 

Fig. 12 - Unitairzation of the scattering amplitude for the qqgi...gn Fock state of the photon by 
the s-channel iteration of the QCD pomeron exchange. 

Fig. 13 - The dominant multipomeron interactions in the deep inelastic scattering. 

Fig. 14 - The multiparticle Reggeons. 
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